Impurity induced low-energy resonances in Bi2Sr2CaCu208+5 
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We study sharp low-energy resonance peaks in the local density of states (LDOS) induced by Zn 
impurities or possible Cu vacancies in superconducting Bi2Sr2CaCu2 0g+5. The measured structure 
of these near-zero-bias resonances is quantitatively reproduced by an extended impurity potential 
without invoking internal impurity states or sophisticated tunneling models. The Zn potential 
extends at least to the nearest-neighbor Cu sites, and the range of order parameter suppression 
extends at least 8 A away from the Zn site. We further show that the local spin susceptibilities near 
Zn impurities increase rather than decrease with decreasing temperature in the superconducting 
state due to the sharp increase of LDOS near the Fermi level. 

PACS numbers: 74.62.Dh,74.25. Jb,74.20.Rp,74.72.Hs 
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Recent scanning tunneling microscopy (STM) 
measurements on the high- Tc superconductor 
Bi2Sr2CaCu208+d- (BSCCO) have revealed the de- 
tailed spectral and spatial structures of sharp resonances 
inside the superconducting gap induced by the impurities 
Zni and Nfi, and by defects such as Cu vacanciesi^ The 
properties of these resonances depend on interactions 
at the atomic scale between impurity potentials and 
the coherent superconducting state. As work progresses 
to probe these defects in the underdoped regime, and 
above T^,, it is essential to understand whether the 
properties in the superconducting state can be under- 
stood within d-wave ECS theory, or whether vestiges 
remain of the exotic correlated state present at higher 
temperatures^iS'S'i Detailed theoretical studies have 
shown that the resonances induced by a Ni impurity in 
BSCCO can be well described by Bogoliubov quasiparti- 
cles scattering off a potential in a d-wave superconductor 
with a small, but discernible, local order parameter 
suppression fS. On the other hand, theories based on 
pointlike potentials have failed to explain quantitatively 
both the spatial structure of the resonances and the 
equally important information yielded by the spatial 
structure of the "coherence peaks", either for Zn 
impurities or Cu vacancies 

Zn impurities induce a sharp resonance at an energy 
close to the Fermi level. Experimental results show that 
these resonances possess a large spectral weight of local 
density of states (LDOS) at the impurity site, whereas 
theories based on pointlike potentials predict that the 
spectral weight at the impurity site is strongly suppressed 
by the unrealistically large strength of the potential. This 
discrepancy in the spatial structure has led to various 
proposals. Some suggest the need for strongly correlated 
models, in which the internal states of the impurity are 
important)^ and others suggest that the spatial struc- 
ture measured by STM is filtered by the surface Bi-0 
layer above the Cu-0 planes. ^'^ However, none of these 
proposals has provided a coherent quantitative account 
for the STM spectra, at both the resonance and the co- 
herence peaks, for both Zn and Ni impurities. In this 
Letter, we show that the effects of Zn impurities over 
the entire spectrum can be quantitatively described by 



spatially extended potentials, and, therefore, render the 
potential model a fully quantitative description for both 
Zn and Ni impurities in BSCCO. A principal, and un- 
expected, conclusion is that the potentials required to 
model the Zn impurity are nonmagnetic and relatively 
weak (< 100 meV), comparable in strength to the po- 
tentials induced by the Ni impurity. Another important 
result is that the order parameter suppression near Zn 
is much more extensive than for Ni. We further report 
that the Cu vacancy^ is well described within such a 
potential model, thus demonstrating that the potential 
model works for all three known defects probed by STM 
in BSCCO. 

The appearance of near-zero-bias resonances is intrin- 
sically important because they greatly influence the bulk 
response to external fields, such as conductivity and spin 
susceptibility. Here we also calculate the local spin sus- 
ceptibilities near a Zn impurity. In the homogeneous bulk 
superconducting state, the spin susceptibility decreases 
with decreasing temperature due to the opening of a gap 
in LDOS around the Fermi level. However, in the pres- 
ence of a sharp near-zero-bias resonance, the local spin 
susceptibilities near the impurity increase with decreas- 
ing temperature, and eventually the values peak at a low 
temperature scale set by the resonance energy. In other 
words, above this temperature scale, the local spin sus- 
ceptibilities inversely scale with temperature, similar to 
the Curie behavior of a paramagnetic ion in an insulator. 
We find such behavior describes the Knight shifts on and 
near Zn in BSCCO. 

We begin by describing the nature of the spatially- 
extended potential that quantitatively describes all of 
the local spectra near the Zn impurity (shown in Fig.^. 
The failure of a pointlike potential to induce the correct 
spatial structure of Zn's near-zero-bias resonance origi- 
nates from the large value of the potential (> 1 eV) re- 
quired to draw states deep into the center of the gap. The 
strong on-site potential diminishes the spectral weight of 
the LDOS at the impurity site, and moves the largest 
spectral weights to the four nearest-neighbor sites. Ex- 
perimentally the largest spectral weights, in order from 
largest to smallest, are at the impurity site, at the second- 
nearest-neighbor sites, and then at the third-nearest- 
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1: A schematic diagram showing the proposed model 
of a Zn impurity on the Cu-0 plane. Parameters extracted 
from fitting to the data are listed on the right. The super- 
conducting order parameters on the dashed links are com- 
pletely suppressed. This model also works for the Cu vacancy 
with Vo = 110 meV, Vi = -50 meV, and 5toi = 19 meV. 
For the Cu vacancy the order parameter was suppressed to 
zero (<5Aij — — A^j) within 8 A of the vacancy, and by 70% 
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-0.7A,w) from 8 A to 12 A from the 
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neighbor sites (See Figs. |21 and 0)1. Such a structure is 
consistent with a non-zero potential at the four nearest- 
neighbor sites. A key feature is that the potential at 
the impurity site is of the opposite sign to those at the 
nearest-neighbor sites, consistent with the structure of 
a screened charge perturbation. The presence of a non- 
magnetic spatially extended potential and a short-ranged 
charge oscillation is consistent also with the results of nu- 
clear quadruple resonance (NQR) measurements on Zn 
impurities in YBa2Cu408 (YBCO)iS These potentials at 
the nearest-neighbor sites are much more effective than 
merely an on-site potential in inducing a sharp resonance 
close to the Fermi level. Unrealistically large values for 
the potential are no longer required. A moderate-range 
local suppression of the order parameter also further 
drives the resonance energy toward the Fermi level. 
Thus a quantitative understanding of the experimental 
data requires careful consideration of the detailed poten- 
tial structure. The same model, but with different values 
for the potentials, including a longer-range order param- 
eter suppression, describes the Cu vacancy. 

We use the Koster-Slater approach to calculate the 
LDOS around impurities AiLiS Homogeneous BSCCO is 
modeled by a one-band tight-binding Hamiltonian, as de- 
scribed in Ref. 0. The Green's function of the homoge- 
neous system is first evaluated with an energy resolution 
5. The Green's function with impurities is then calcu- 
lated by solving the Gorkov (Dyson) equation. The im- 
purity potential consists of site-diagonal matrix elements 
at the impurity site (Vq) and at the nearest-neighbor sites 
(Vi), and off-diagonal matrix elements corresponding to 
modifications of the hopping matrix elements {5t) , and of 
the superconducting order parameters (^A) as described 
in Ref. (3 for the Ni impurity. The full Hamiltonian that 
describes the one impurity problem takes the following 
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where i and j label the lattice sites (the impurity re- 
sides at site 0), and cr labels spin. The hopping ma- 
trix elements, Uj = {148.8,-40.9,13,14,-12.8} meV, 
are based on a one-band parameterization of the angle- 
resolved photoemission dataJi^ The superconducting or- 
der parameters, A^j, of the homogeneous system are only 
non-zero on the bonds connecting two nearest-neighbor 
sites, Ai^j+a; -Aj^j+y = Ao/4, where Aq is the gap 
maximum. The momentum-dependent order parameter 
resulting from these A^^ has d-wave symmetry, Ak — 
(Ao/2)(cos kxL — cos kyL), where L is the lattice spacing 
between two Cu atoms. 

To construct the effective potential for the Zn impu- 
rity, we relate important features in the LDOS data to 
different parts of the potential structure. First of all, it 
is known from STM that there is microscopic inhomo- 
geneity of the LDOS gap, which appears to be related to 
the local doping concentration.^^ Judging by the exper- 
imental LDOS spectrum at a site within the same local 
patch as the Zn impurity, but reasonably far away from 
the impurity, we determine that the gap maximum Aq is 
about 32 meV in Ref.Q. The local chemical potential /i is 
shifted with respect to the optimal doping value (—130.5 
meV)i^ by —20 meV, which sets the Van Hove peak to be 
about 14 meV below the Fermi level. The energy resolu- 
tion 5 is about 2 meV, roughly judging by the widths of 
the coherence peaks and of the Zn resonance. The range 
of order parameter suppression is then determined by the 
closing in of the local gap edge. Because the amplitude of 
the remnant coherence peak is strongly suppressed and 
the appearance of mid-gap resonances contributes addi- 
tional spectral weight inside the gap, our fit for the gap 
edge is not as quantitative as for other features of the 
resonance. Nevertheless, the minimum range of the or- 
der parameter suppression is determined to be about 2 
lattice constants away from the impurity. That is, the 
order parameter is completely suppressed within a circle 
with a radius of approximately 8 A. This is a much more 
extensive order parameter suppression than seen around 
the Ni impurity,^ and suggests that Zn is more destruc- 
tive to local superconductivity. 

The matrix elements Vb, Vi and J^oi are then deter- 
mined based on the resonance energy and the peak am- 
plitudes at three different sites (the impurity site, the 
1st and the 2nd nearest-neighbor sites). Although all 3 
parameters jointly determine the spatial structure, Vi is 
the most dominant component. Vq is used to obtain the 
correct resonance energy and St^i is used to adjust the 
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FIG. 2: LDOS spectra per unit cell at various sites near a 
Zn impurity. The left panels show the spectra in the "nor- 
mal state" (simply setting Aij = 0) using the same impurity 
parameters. The right panels show the spectra in the su- 
perconducting state. Solid lines show the calculated results, 
(a) and (b) are the spectra at a site far away from the Zn 
impurity, (c) and (d) are the spectra right at the impurity 
site, (e) and (f) are the spectra at the 1st nearest-neighbor 
sites, (g) and (h) are the spectra at the 2nd nearest-neighbor 
sites. Open circles (o) show the STM differential conductance 
datai- The data was rescaled by a constant factor identical 
for all the spectra. Note that the effective potential in the 
"normal state" could be very different. 



ratio between the peaks. Our procedure provides simi- 
lar spatial structures if the sign of both Vq and Vi are 
reversed - we can only determine that they have oppo- 
site sign relative to each other li^ Motivated by the NQR 
experiments that suggest electron charge accumulation 
on the nearest-neighbor sites of Zn in YBCO, we choose 
the potential energy on the nearest-neighbor sites to have 
negative values (attractive to electrons). 

Figure |2] shows the detailed comparisons of our calcu- 
lations with data. In general, we have a good fit for the 
spectra in a wide energy range inside the gap. The fit is 
less satisfactory for the energy range beyond the coher- 
ence peaks. Therefore, the determination of the gap edge 
is not as precise as the determination of other parameters 
based on the resonance energy and peak spectral weights. 
We note that for every midgap resonance one expects a 
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FIG. 3; (color) Spatial structure of the LDOS near a Zn 
impurity at —2 meV in (a) 2D and (b) 3D view. We have 
spatially distributed the LDOS at each unit cell according to 
a normalized Gaussian with a width of half the lattice spac- 
ing (L). The LDOS is shown in logarithmic scale. The Cu-O 
lattice is rotated by 45° so that the horizontal axis is aligned 
with the superlattice modulation. To demonstrate the effect 
of the band structure, we carried out similar calculations for 
tight-binding bands with (c) only the nearest-neighbor hop- 
ping ({ti, At, Vo, Vi} = {150, 10, 100, -100} meV) and with (d) 
only the nearest- neighbor and the second-nearest-neighbor 
hopping {{ti,t2,fj., Vo, Vi} = {150, -40, -150, 75, -75} meV). 
The chemical potential is chosen so that the Van Hove singu- 
larity is fixed at the same energy. The impurity potential is 
also slightly altered to keep the resonance energy at the same 
place. 



feature on the nearest-neighbor sites which is on the op- 
posite side from that on the impurity site itselfAi^ For 
these spatially-extended potentials, however, the spec- 
tral weight of these resonance peaks is highly reduced 
and likely not visible. 

The spatial structure of the resonance peak, suitable 
for direct comparison with Ref. 1, is shown in Fig.|3| Note 
that the spectral weight is most peaked along the (1,1) 
direction of the Cu-O lattice in the near region around 
the impurity. Far from the impurity, the tails align with 
the Cu-O bond direction (This is also the case for Ni, 
but the signal is much weaker.) As our fit of the local 
potentials was constrained entirely by features near the 
impurity, obtaining good agreement with measured data 
for this orientation of the resonance tail suggests that 
we have a reasonable description of the BSCCO band 
structure near the Fermi level. For example, in a sim- 
ple nearest-neighbor tight-binding model we find that the 
long-ranged tail of the resonance incorrectly follows the 
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FIG. 4: The local spin susceptibilities at the Zn site (solid 
line), at the 1st nearest-neighbor sites (dot-dashed line), at 
the 2nd nearest-neighbor sites (double dot-dashed line), and 
at a faraway site (dotted line). The susceptibilities near Zn 
show a Curie-like behavior in an intermediate temperature 
range (above 10 K and below the superconducting transition 
temperature). The inset shows the corresponding susceptibil- 
ities near a Ni impurity, which simply increase with temper- 
ature. 

diagonal of the Cu-0 lattice [Fig.|3Jc)]. The tail begins to 
rotate to follow correctly the Cu-0 bond direction when 
the second- nearest- neighbor hopping is added [Fig.lS^d)], 
and the proper orientation is obtained in the full tight- 
binding model for BSCCO. We have verified that this 
defect model, with the different potentials reported in 
Fig. n also describes well the more limited Cu vacancy 
data in Ref. 0- The order parameter suppression even 
more extended around the Cu vacancy than around the 
Zn impurity; our model parameters are listed in the Fig.^ 
caption. 

To show that the bulk response to fields can be dra- 
matically changed by the appearance of the near-zero- 
bias resonance, we calculate the local spin susceptibilities 
(suitable for determining the Knight shifts) around the 
impurity using our results for the LDOS. The local spin 
susceptibility at site j is defined as 

X. = ^ / <h;Mu;)fic.)[l^fiu;)], (2) 
where /xs is the Bohr magneton, ks is the Boltzmann 



constant, T is temperature, Aj{uj) is the LDOS at site 
j, and /(w) is the Fermi-Dirac distribution function. We 
assume g — 2. The results for both Zn and Ni impurities 
are shown in Fig. 01 For Zn the local spin susceptibil- 
ity increases as the temperature decreases, and there is 
a maximum at a very low temperature, which is deter- 
mined by the resonance energy. The strongest response 
is at the impurity site, and the signal strengths at the 1st 
and 2nd nearest-neighbor sites are approximately equal. 
Ni, which has a much higher resonance energy, does not 
show this low-temperature maximum, and thus does not 
appear to be paramagnetic in this type of experiment. 
Thus we find that the resonance induced by a completely 
nonmagnetic impurity potential for Zn generates a signal 
in spin susceptibility experiments that closely mimics the 
expected signal from a free spin (as suggested in Ref.0). 
This result may have implications for measurements of 
spin susceptibilities in YBa2Ca307_5fi£ 

We have provided a reasonable nonmagnetic potential 
model that quantitatively reproduces the spatial struc- 
ture of the Zn impurity resonance and the spatial struc- 
ture of the coherence peaks near Zn. We have verified 
that the potential model also produces accurate results 
for the Cu vacancy. The model is based on the same 
BSCCO electronic structure as used previously to de- 
scribe Ni, thus demonstrating that a consistent picture 
is possible that quantitatively describes all three defects 
in BSCCO. Surprisingly small effective potentials are re- 
quired to describe the Zn resonance, even though it oc- 
curs near zero energy. The order parameter suppression 
we find near Zn is more extensive than near Ni, but not as 
extensive as near the Cu vacancy. We have also demon- 
strated that the resonance near Zn would respond in spin 
susceptibility measurements in a way that would mimic 
a free spin, whereas Ni would not. Finally, the potential 
for the vacancy, and to a lesser extent Zn, scatters quasi- 
particles most effectively for momentum transfer (tTjTt). 
This may broaden quasiparticle signatures significantly 
at the (tt, 0) points of the Brillouin zone. These results 
provide a foundation for better understanding of mea- 
surements that are likely to be made on these defects in 
the more exotic "pseudogap" state of BSCCO at higher 
temperatures. 
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